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Abstract 

The irrotational motion of a compressible inviscid fluid is studied in the field of analogue 
gravity, where its metric is compared to that in general relativity [1] , a fluid analogue of an evap- 
orating black hole [2] has been realized experimentally [3], and there are symmetries related to 
the standard model [3]. Here we show the analogy also extends quantitatively to electromag- 
netic and quantum mechanical phenomena. We discuss a candidate model to account for the 
number and precision of these analogies. 



1 Introduction 

We will show that Euler's equation for a com- 
pressible fluid has irrotational solutions which 
superficially resemble smoke rings. Their equa- 
tions of motion through the fluid are completely 
classical, and the resulting trajectories obey sta- 
tistical equations which are identical to those of 
the Copenhagen interpretation of quantum me- 
chanics. A similar mechanism accounts for the 
emergence of quantum phenomena from classical 
motion in experiments where a droplet is made 
to bounce on a liquid surface, exhibiting quan- 
tised energy levels, single-slit diffraction, double- 
slit diffraction and tunnelling [51 El El [HI IS] - Some 
of the ring-like solutions are chiral. Opposite chi- 
ralities attract and like chiralities repel with an 
inverse square force which has the same fluid dy- 
namic origin as the effect used for degassing oil 
by subjection to ultrasonic vibration [TDJ §4.4]. 
The interaction obeys Maxwell's equations and 
its strength is characterized by a fine structure 
constant, a < 1/49. 

Consider a compressible inviscid fluid such as 
the air in the idealization it is continuous and 



has no viscosity or thermal conductivity. Such a 
fluid obeys Euler's equation [10] 

d i + (uV)u = " l P VP (1) 

where the pressure P(x, t) is a function of the 
density p(x, t) and dp/dt = — V(pu). At low 
amplitude this reduces to the wave equation 

where c 2 = |p. In addition to describing sound 
waves, this equation has cylindrical solutions, 
the simplest of which is present in the surface 
waves just after a raindrop has hit a puddle. In 
cylindrical coordinates (r, 9, s) 

i = A e -i("ot+me-k sS ) j m ( fcrr ) ( 3 ) 

where y = 1 + -Re(£) and p Q is the mean density, 
J m is a cylindrical Bessel function of the first 
kind, m an integer and u 2 = c 2 (k 2 + k 2 ). It may 
be verified by substitution that this obeys the 
wave equation. The flow pattern is irrotational, 
that is, § u.dl = for any closed path, as shown 
in the appendix. 
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Figure 1: Sketch of an eddy, described by (|3| 
with m = 1 and k s = 0. The compressions 
(red, schematic and not to scale) and rarefac- 
tions (green) rotate as shown. Note the contrary 
direction of fluid flow near the rarefaction. 

Figure [l] sketches the m — 1 solution. There 
are families of solutions, superficially resembling 
smoke rings, which can be obtained by curving 
it into a torus. Two such solutions, which we 
call 'sonons', are sketched in Figure [2j 




Figure 2: Two sonons obtained by curving the 
eddy of Figure [I] into a torus. The compressions 
(red) and rarefactions (green) rotate as shown, 
(a) i?io (see equation |6| (b) Rn, obtained by 
adding a twist. This solution is chiral and the 
motion resolves into a spin about the axis. 

In the case m = or 1, the Appendix shows 
that the associated density pattern can be writ- 
ten to a good approximation as a sum of spher- 
ical harmonic solutions to the wave equation in 
the coordinates of figure [3}- 

f = Mt) ^mn(x) (4) 

where 

vjj o = Ae~ iuJot (5) 

R mn = 1^ e-^ me '- n ^ ]m {k r a)k r R d<P (6) 
Jo 

where n is an integer and j m is a spherical Bessel 
function of the first kind. 




Figure 3: The coordinates used in ([6]) 

1.1 Dimensions of a sonon 

If the eddy in Figure [T] is terminated on two pla- 
nar surfaces, they will attract one another due 
to the reduced Bernoulli pressure. This must 
be due to an attraction along the length of the 
eddy, similar to the attraction in a vortex. 

When an ordinary vortex is curved into a 
smoke ring, this force is balanced by Mag- 
nus forces (like the lift of an aircraft wing) 
as the structure moves forward through the 
fluid [10J. However a sonon cannot experience 
Magnus forces because it is irrotational, and con- 
sequently its radius will shrink, causing the am- 
plitude A in (|5| to grow due to the conserva- 
tion of fluid energy. Nonlinear effects will halt 
the shrinking before A reaches about 1 since the 
density cannot become negative. 




Figure 4: (a) The fluid flows (blue) near a Rn 
sonon reinforce in the y direction at the instant 
shown, (b) The nonlinear region of reinforce- 
ment approximates to a cylinder. If it has di- 
ameter is aR a , length 2bR Q and mean flow speed 
| (3c, a reasonable estimate isa~6~/3~l 

Very near a Rn sonon, shown in Figure |4| the 
main component of the flow derives from the ra- 
dial flow speed due to an eddy with m = 1, 
which is u r oc ^J\{k r r) (see the appendix). 
These flows reinforce provided the diameter of 
the sonon is less than the first zero of u r , which 
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is 2 x 0.92. Neglecting curvature and nonlinear- 
ities, the lowest energy (unexcited) state has, as 
an upper estimate, 



R 



< 



0.92 



(7) 



1.2 Spin symmetry 

The integral in ^ can be factorised into R mn = 
[£ r .$ r (r)][£ e .$ fl (0)][B .$ (0)], where B, and 
are vectors. 

At large distance, the components of $ r can 
be obtained by expanding the Bessel function at 
large r and neglecting small terms 



$ r = - [ sin( k r r), cos( k r r) 
r 



Notice that 







where we have defined $ = $ r $e$0. The other 
eigenvector, [sin(Avr), — cos(fc r r)], has an eigen- 
value of the opposite sign. 

By inspection of ^ there is an eigenvector 
®4> = [e~ i ™ t> ,e in<t> ], which obeys 



1 \ <9$ 
/ I Q _ 1 I — = n$ 



Reversing the sign of n corresponds to revers- 
ing the spin direction, and it produces the other 
eigenvector with the opposite sign of eigenvalue. 

In the 9 direction, the contributions from the 
near and far sides of the sonon produce wave 
patterns which propagate in opposite directions 
at the calculation position in Figure [3j The two 
orthogonal components are = [$ + , $~] where 



-imd 



±e ime . Thus 



1 

1 



<9$ 



Notice that these matrices are the same as the 
Pauli spin matrices. The appendix discusses a 
possible physical interpretation of the associated 
spin-half symmetry, and also the analogy with 
the Dirac equation. 



1.3 Boundary condition 

Taking one of the components of $ r in Q, the 
kinetic energy | f pu 2 can be calculated in the 
usual way in fluid dynamics. The value must 
be doubled to account for the potential energy 
of motion. The total energy inside a sphere of 
radius R is (see the appendix) 



E 



fluid 



tiR 2 R p a c 2 



(9) 



This approximates to the energy in the fluid 
motion of sonon since A and B r i are of order 1 
and the others diminish rapidly with radius. 

Even without container walls, there will still 
be a boundary condition due to other sonons, 
since the energy in (|9| will be reduced if they 
align themselves so they interfere destructively 
at large distance. The coupling mechanism 
which allows the lowest energy state to be 
reached was first reported in pendulum clocks 
by Huygens in 1665 [TT] - see the appendix. 

Figure [5] shows two Ru sonons which are 
aligned so they interfere destructively at large 
distance. The spherical Bessel function j\{r) 
changes phase by ir compared to a freely propa- 
gating wave when going through the centre, and 
it follows that destructive interference at large 
distance means there must be constructive inter- 
ference in the region between the sonons. There 
is consequently an antinode of density on the 
mirror line. Notice that the spin vectors in Fig- 
ure [5] are opposed. 



(a) ffp 



(b) ^ 




M 



M 




Figure 5: Pairs of Ru sonons of (a) opposite 
(b) the same chirality, aligned so their waves in- 
terfere destructively at large distance on an ex- 
tended line joining them. There is an antinode 
of density on the mirror line M. 
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Suppose the sonons in a region are brought 
closer together to form a denser body where the 
typical distance R between them is reduced. If 
the process is adiabatic then, from (|9]), RR 2 will 
remain constant, and consequently the sonons 
will become enlarged, with a greater value of R Q . 
These enlarged structures are valid solutions to 
Euler's equation provided their frequencies are 
also redshifted, since if £(x, t) is a solution then 
it may be verified by substitution into (|TJ that 
((ax, at) is also a solution where a is an arbitrary 
scale factor (noting that the velocity u is unaf- 
fected by the scaling). A similar redshift near a 
massive body is studied in analogue gravity, as 
discussed further in the appendix. 

2 Lorentz covariance 

From ([9]), most of the fluid energy of a sonon 
is at large distance, where the motion obeys the 
wave equation because it is low amplitude. Since 
the wave equation is Lorentz covariant, so is this 
motion. Thus, if £(x, t) is a solution then £(x', t') 
is also a solution where x' = j(x — vt), y' = y, 
z' = z, t' = 7(t — vx/c 2 ) and - = \A — v 2 /c 2 . 

Very near the sonon, the amplitude of mo- 
tion is not negligible and the Lorentz symmetry 
might, in principle, be perturbed by the non- 
linear term (u.V)u in Q. However, consider 
adding a constant velocity to all motion. Re- 
placing u by u + v where v is constant gives 
a perturbation e = (v.V)u. This perturbation 
vanishes over a cycle because the fluid motion 
is oscillatory and returns to its starting-point, 
J u dt = and so 



e dt 



V.V)u dt 







Expectation values are usually defined to con- 
verge on the long term measurement. For preci- 
sion we define an 'ordinary' expectation value 
where this measurement encompasses one or 
more complete cycles. It follows that no devia- 
tions from Lorentz covariance can be discovered 
by measuring any ordinary expectation value. 

Another potential perturbation arises because 
the speed c = a/ dPf dp varies over the cycle. 



However, in the absence of external influences, 
the speed of a propagating wave is unperturbed 
because its leading edge cannot travel faster 
than the low amplitude speed, as is known for 
sound waves which obey similar equations [TU] . 

Finally, a sonon might be perturbed by the 
mean variation in c 2 = |p. However, the effect 
is third-order, as can be seen by expanding P(p) 
in terms of s — (p — p )/p : the mean of ^ 
is only perturbed by the odd-order terms, start- 
ing with s 3 . Even if the effect were to rise above 
noise levels, we will see it would be manifested as 
an adjustment to the effective mass of the quasi- 
particle, which would be difficult to distinguish 
from other effects such as thermal energy. 

It follows that all ordinary expectation values 
are Lorentz covariant to a very good approxima- 
tion at all amplitudes. 



3 Equations of motion 

From Q, V 2 ip Q = and ^ Q = -ufyj) . We 
have seen the sonon is Lorentz covariant to a 
very good approximation, so this must be a 
special case of a Lorentz covariant relationship 
which is velocity-independent, namely 



d 2 i) 



c 2 V 2 ^ 



(10) 



which is analogous to the Klein-Gordon equation 
for a relativistic particle. At low velocity, (10) 



can be approximated by defining ifj = e luJot ifj' 



and neglecting JW 



giving 



vy (11) 



where, in a change of units to match convention, 
we have defined m = % and V is a constant. 
This is analogous to the Schrodinger equation 
for motion in a constant potential. 

To complete the description, applying a 
Lorentz boost to the density pattern R mn in ^ 
and neglecting its time dependence (which is sec- 
ond order in v) gives the same equation with r 
replaced by r' where r' x = j(r x — vt), r' y = r y 
and r' z = r z . We will refer to this as x(x). Its 
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peaks and troughs also move through space at 
the velocity v, and consequently the sonon will 
remain aligned with them if it is unperturbed. 
This coherence is likely to be reinforced by non- 
linearities in the medium. We will first examine 
the coherent case, and then examine the effects 
of decoherence. 

An expression for the velocity v in terms of 
V> follows by applying a Lorentz boost to (|5p, 
giving ip = Ae where S = k.x — cut and, from 
the Lorentz transformation, u = 7Co> and k = 
(7w /c 2 )v. These combine to Kk = mv where 
we have defined m = ^ . This leads immediately 
to 

y=ll m (^) (12) 

The significance of this expression can be under- 
stood by following Bohm [12], who rearranged 
the Shrodinger equation into the form 



+ V 







(13) 



where v is given by ( 12 ) . Bohm noticed that ( 13 ) 



is the standard equation for a conserved quantity 
moving at velocity v, and concluded that \ip\ 2 
must be the probability (averaged over nearby 
trajectories) that a hypothetical particle mov- 
ing at velocity v will pass through a given posi- 
tion at a given time. These particle trajectories 
are the same as those of sonons. However the re- 
sulting de Broglie-Bohm or 'pilot wave' model of 
quantum mechanics [T2"j IT5] describes only part 
of the motion of sonons because it omits the car- 
rier waves. 

Suppose, conterfactually, that a sonon 
changes state and becomes de-localized in some 
unspecified way, and then suddenly rematerial- 
izes just before any measurement, without af- 
fecting the trajectory. It follows from the above 
that the probability of rematerialization at any 
given position must be \ip\ 2 . This is analogous to 
the Copenhagen interpretation of quantum me- 
chanics, where the wavefunction and its collapse 
obey identical equations. 

The foregoing can be described in the termi- 
nology of telecommunications theory. At veloci- 
ties significantly smaller than c, the wavelength 



of x (the carrier wave) is significantly shorter 
than that of ip (the modulation). Note that ip 
must be complex-valued to describe both ampli- 
tude and phase. 

A similar emergence of quantum phenomena 
from completely classical motion is exhibited by 
small droplets which are made to bounce on a 
bath of the same liquid by oscillating the con- 
tainer vertically [5J. The droplet's horizontal 
motion is determined by nonlinear interactions 
with the surface waves, which obey the wave 
equation and £1X6 cL two-dimensional analogue of 
the disturbance £ = x4>- The principal differ- 
ences are that the effective speed c might be 
modified near the bouncing droplet due to the 
mass of the droplet and amplitude-related ef- 
fects, and the experiment is a driven dissipative 
system rather than a lossless one. The videos 
show the droplets maintaining phase coherence 
with the surrounding waves due to nonlineari- 
ties. Single-slit diffraction, double-slit diffrac- 
tion, unpredictable tunnelling and quantized en- 
ergy levels are observed (HI |7J El [9] . 

A sonon may be modelled £IS db classical os- 
cillator of high Q. In interactions between such 
oscillators, Mead shows that the sum of the fre- 
quencies is conserved, since otherwise the inter- 
action would average out to zero [H]. It follows 
from the definition E = hu (note, E is not the 
same as energy in the fluid motion Efi uid referred 
to above) that E is conserved. The conservation 
of momentum, defined by p = hk., follows from 
the conservation of E and Lorentz covariance. 

It follows from the above that the equations of 
motion for a sonon are completely classical, but 
they reduce to the established equations for a 
relativistic quantum mechanical particle of mass 



if the carrier wave is neglected. 



4 Long range force 

In the bouncing droplet experiments [Zj, the 
edge of the container can be modelled as 
an image droplet bouncing antiphase, leaving 
the meniscus undisturbed. Stroboscopic pho- 
tographs show a droplet moves in a straight line 
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until it approaches the edge, when it is deflected 
away. We interpret this to indicate a repulsive 
force from the antiphase image. Conversely, in- 
phase droplets form crystalline lattices, indicat- 
ing an attractive force. 

A similar phenomenon in three dimensions is 
used for degassing oils by subjection to ultra- 
sonic vibration. As the bubbles expand and con- 
tract in response to the pressure waves, they in- 
duce flows in the fluid. When two bubbles are in- 
phase, no fluid will cross the mirror line between 
them, from symmetry, but the flows parallel to it 
reinforce, resulting in a reduced Bernoulli pres- 
sure and a force of attraction. The bubbles con- 
sequently merge and rise to the surface. 

The calculation of the force, adapted from 
Faber's textbook [10j §4.4], is as follows. If a 
vacuum cleaner hose operates in the wind, it 
will ingest momentum along with the air par- 
ticles, and experience a force F = p UQ\ where 
Qi is the flow and U the wind velocity. If the 
wind is due to a similar hose with flow Q2, then 
U = Q2I '(4-7rr 2 ), giving an inverse square force 
of attraction 



- Po 



Q1Q2 



(14) 



The direction of flow, and of the force on both 
hoses, will be reversed if one of them is set to 
blow (with a baffle for spherical symmetry). If 
both are set to blow they will attract again. 
More generally, oscillatory motion results in an 
attractive force if it is in-phase, and a repulsion 
if it is antiphase. The phase alignment is most 
easily assessed from the flow across the mirror 
line between the sources. 

Turning to Ru sonons, their flows are shown 
in Figure [4] and their preferred alignments in Fig- 
ure [5] From the mirror symmetry in[5|a) there is 
no flow across the mirror line and consequently 
sonons of opposite chiralities will attract one an- 
other. In [5^b) the direction of flow is reversed 
and like chiralites repel. The magnitude of the 
force is calculated in the appendix by substitut- 
ing the idealised fluid flows of Figure |4|(b) into 



(a calculation which is only approximate) 

F " w c sin2(C) % (15) 

where a, b and /3 are defined in Figure [4j and 
C = ^k 2 r E? + 1. Substituting the estimates a « 
b « « 1 and k r R < 0.92 from Q gives 

1 he , , 

F < (16) 

~ 49 r 2 1 ' 

This inverse square force among Ru sonons and 
their chiral twins is analogous to the Coulomb 
force among electrons and positrons. The fine 
structure constant a < 1/49 matches the elec- 
tromagnetic value to the accuracy of calculation. 

4.1 Radiation 



The interaction in ( 16 ) can be generalized by ap- 



plying a Lorentz boost. We can borrow from the 
usual extension of the Coulomb force by requir- 
ing it to be Lorentz covariant, which produces 
Maxwell's equations, as shown in the appendix. 

Maxwell's equations have solutions corre- 
sponding to propagating waves. We now show 
these waves are emitted when a sonon is acceler- 
ated. Without loss of generality we will examine 
one component of <3? in ([8]). 

If a sonon is made to oscillate so its position 
is x = x Q sin(f2t), where x D is small, then substi- 
tuting into Q gives (see the appendix) 



where 



Xv 



AR n 



B 1 e 



k r oc q 



sin 



n (t 



(17) 



and B is the angle from the direction of oscilla- 
tion. Here ip w describes amplitude modulation 
of a carrier wave Xw, but in general there will be 
phase modulation as well for larger amplitude 
motion. 

The momentum of the wave can be calculated 



(14) and calculating the cylinder's acceleration 



by writing (|17|) as a sum of propagating waves 

-iu> (t+r/c) 



+ e 



(18) 
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These propagating solutions can be extended 
all the way to the particle (more precisely using 
a sum of Hankel functions). Extended in this 
way, they comprise the entire solution. Their 
sum must have the same energy and momentum 
as the particle, at all velocities, which enables 
a relationship between the energy E w and mo- 
mentum p w of the waves to be calculated 



E 



Pu 



w 
C 



See the appendix for the calculation. This also 
applies to the radiating waves described above, 
which are a superposition of such waves. Notice 
the analogue with light, which carries the same 
momentum. 



4.2 Constructive interference 

In very shallow water on a sandy beach, the 
wave speed increases with depth, and a wave 
travelling close behind another will propagate in 
the deeper water of the wave in front and travel 
faster, catching it up [TO] . 

A similar effect will occur in the waves un- 
der consideration due to the speed c = a/ dPf dp 
varying over the cycle. The trailing wave will be 
perturbed by the wave in front, which remains 
in phase over a long time, and therefore it does 
not cancel over a cycle. This effect is not an or- 
dinary expectation value (as defined in Q and 
therefore it is not precisely Lorentz covariant. 

There is an analogy with light waves, which 
also tend to come into coherence, such as in 
lasers. The underlying components of the waves 



in (18) will interfere constructively if one over- 



takes the other by about half a wavelength of the 
carrier wave. For comparison, an oscillating elec- 
tron would have about 10 12 half-wavelengths per 
metre in its carrier wave. On a crude model, two 
waves of modulation will reach coherence after 
about lm if the speed difference is Ac ~ lCT 12 c, 
which corresponds to about 0.0003ms" 1 . 



5 Causality and localization 

Euler's equation is deterministic, and therefore 
chance and dice cannot be involved at a fun- 
damental level in the behaviour of sonons in a 
closed system. The difficulty of predicting out- 
comes arises, among other causes, from errors or 
uncertainties in the initial measurements, which 
are magnified in processes such as diffraction and 
tunnelling. 

In a closed deterministic system, its state at 
some time t completely dictates its state at 
some other time t. This is true even when t > t, 
which could be interpreted non-causally. Such 
solutions are valid because Euler's equation is 
symmetric under time reversal. 

In practice, these coherent solutions are easily 
disrupted by influences from outside the closed 
system. In particular, the coherence between a 
sonon and its surrounding waves is fragile. The 
process of decoherence is most easily pictured 
in the droplet experiments, where a disturbance 
causes a droplet to bounce out of its local wave 
trough, magnifying the effect of the perturba- 
tion. 

Nevertheless, the time-reversed motions are 
valid solutions if decoherence can be controlled. 
See Mead's analysis [U] from first principles 
on the interactions between classical oscillators 
through forward and time-reversed waves which 
are solutions to the (time-symmetric) wave equa- 
tion. Interacting sonons obey similar equations. 
Provided phase coherence can be maintained 
during the transition, Mead shows that the in- 
teraction produces effects which are closely anal- 
ogous to the quantisation of the photon, among 
other observed phenomena. 

Mead's model is related to Cramer's more gen- 
eral transactional interpretation of quantum me- 
chanics [13], which was designed to be consistent 
with experiments on Bell's inequality [TBI, [T71 [18] 
by exploiting the property that the equations are 
symmetric under time reversal. Cramer's model 
is completely local and time reversal symmetric, 
symmetries which are shared by Euler's equa- 
tion. It might be possible to interpret the mo- 
tion of sonons (still consistent with experiments 
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on Bell's inequality) in a different way, without 
needing to invoke the time-reversal symmetry, 
by exploiting the fact that spin-related informa- 
tion is carried by the carrier waves. This trans- 
mission of information is not usually considered 
in the interpretation of the experiments. See 
fTOl for further discussion. 



6 Conclusion 

The irrotational motion of a compressible invis- 
cid fluid is studied in the field of analogue grav- 
ity, where there are phenomena related to gen- 
eral relativity and the standard model [TJ [2J, [3J H] . 
The main contribution of this paper has been 
to extend the analogy to quantum and electro- 
magnetic phenomena. We introduced the sonon, 
a quasiparticle with a twist, which behaves like 
a relativistic quantum mechanical particle with 
spin- 1 symmetry, and which experiences an in- 
teraction, due ultimately to Bernoulli forces, 
that obeys Maxwell's equations. 

The number and precision of these analo- 
gies suggests there might be a wider explana- 
tion. One candidate arises from an argument 
advanced by Einstein in 1920 [2D]. Unlike linear 
motion, acceleration and rotation are intrinsi- 
cally discoverable. For example, a gyroscope's 
axis tends to remain aligned with the distant 
stars. Believing all interactions to be local, Ein- 
stein reasoned these correlations between distant 
objects must be mediated by a substance occu- 
pying the space between them, which he called 
a "medium for the effects of inertia" or "ether of 
the general theory of relativity". By requiring 
the medium to be consistent with special rela- 
tivity he deduced that it cannot have any me- 
chanical properties if it is a solid or a semi-solid. 

A fluid medium was not considered at the 
time because it cannot support transverse waves, 
which were thought to be needed to model po- 
larised light. However, this assumption cannot 
be justified because a compressible inviscid fluid 
supports, not transverse waves, but waves of 
modulation which obey Maxwell's equations. 

These observations are consistent with the hy- 



pothesis, which we now advance, that Einstein's 
inertial medium behaves as a nonrelativistic 
barotropically compressible inviscid fluid. The 
hypothesis is testable because it predicts all ob- 
servable quantities are completely determined 
by the solutions to Euler's equation for such a 
medium. 

The observable - and observed - quantities 
described above in general relativity, particle 
symmetries, quantum mechanics and electro- 
magnetism can be derived from this single hy- 
pothesis. In addition, it predicts that all ordi- 
nary expectation values (as defined in Q must 
be Lorentz covariant to a very good approxi- 
mation at all amplitudes. This is consistent 
with the original paper on special relativity [2TJ , 
which specifies that the postulates apply to the 
speed of light in empty space and the phenom- 
ena of electrodynamics and mechanics (we inter- 
pret this to refer to the corresponding ordinary 
expectation values). However, our hypothesis is 
not consistent with a possible extension to the 
original postulates, in which all motion is pre- 
sumed to be Lorentz covariant. In particular, 
there are rare expectation values which do not 
average over a cycle but still have observable 
consequences. These are not precisely Lorentz 
covariant and, as described in 94.21 such devia- 



tions contribute to, or may be responsible for, 
the coherence effects observed in lasers. 

The calculation of the fine structure constant 
might be improved by computer simulation and 
compared to the observed electromagnetic value. 
The families of quasiparticles and the redshift 



near a massive body in £1.3 might be studied 
further and related to analogue gravity. 

I warmly thank Andrew McLachlan, Keith 
Moffatt, Graziano Brady, Peter Landrock, and 
especially Robin Ball and Ross Anderson for im- 
mensely helpful comment and discussion. 
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Appendices 

A Symbols 



Symbol Meaning 



(x,t) 
(r,6,s) 



P(x,t) 
u(x,t) 
p(x,t) 

Po 

£(x,t) 



jm(r) 

Rmn 
V 

7 

^(x,t) 
(w,k) 



Cartesian coordinates 
Cylindrical coordinates 
Angle around torus 

Fluid pressure 

Fluid velocity 

Fluid density 

Mean fluid density 

Excess density 

denned by pj p Q = 1 + i?e(£) 

Speed of propagation 

denned by c 2 = dPj dp 

Cylindrical Bessel function 

of the first kind 

Spherical Bessel function 

of the first kind 

Spatial dependence of a sonon 

Velocity of Lorentz boost 
Lorentz factor 1 / yT— v 2 /c 2 



Copenhagen wavefunction 

Carrier wave £ = 

Angular frequency, wavevector 



B Flow pattern of an eddy 

The article discusses a solution to Euler's equa- 
tion at low amplitude, given by 

where here, and throughout this appendix, the 
symbols have the meanings in the table above. 

At low amplitude the quadratic term in 
Euler's equation can be neglected, giving 
p du/dt = —VP. Substituting the definition 
c 2 = dP/dp gives p D du/dt = c 2 Vp or 



u 



c 

Po 



Vp dt 



(19) 



This flow is irrotational [TO], as can be seen 
by rearranging ( 19 ) into the form u = V0 where 
= c 2 J s dt where s = (p 



i 



Po)/po- 



Cross-section 



I (• { •)• 1 



• Centre of eddy 

• Stagnation point 

• Fluid flow 

^ Density maximum (+) 
or minimum (-) 



Figure 6: Flow pattern of an eddy (equation [3] 
where m — 1 and k s = 0) at the instant t — 

Figure [6] sketches the flow pattern near such 
an eddy with m = 1, which can be calculated by 



expressing (19) in cylindrical coordinates 



c 2 dp 
Po dr 

c 2 dp 



dt 



dt 



p Q r 36 

Substituting into ^ and integrating gives 

c 2 d 

A — -7^J m (k r r) sin(a; t + m6 — k s s) 
<jj or 



C 2 Til 

ug = A — — J m (k r r) cos(u t + m6 — k s s) 
u r 

In Figure [6] the x axis is aligned with 6 = 
and the flow pattern is sketched at the instant 
t = s = 0. The flow is stationary at the 
points marked with a small dot on the diagram. 
On the x-axis this stationary flow occurs when 
Ji(k r r) = 0, or k r x = {3.83,7.01,10.71..}. On 
the y-axis it occurs when (d/dr)J\(k r r) = 0, or 
k r y = {1.84,5.33,8.53..}. 
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C A sonon as a sum of 
spherical harmonics 

A sonon is obtained by curving an eddy ^ into 
a ring. This is a good description at small r, but 
the curvature of the s axis makes the analysis 
more difficult at larger r. We now describe a 
solution to the wave equation which behaves in 
the same way at large distance, and uses more 
tractable coordinates. 




Figure 7: The coordinates used in (20) 



In the coordinates of Figure [7], we begin by 
showing that 



J m (r) 



/oo 
j m (cr) ds (r -> oo) (20) 



Expanding the spherical Bessel function at 
large distance, the dominant term is 



1 / 7T 

— sin a — m— 
a V 2 



r — > oo 



The integral in ( 20 ) can be evaluated by Fres- 



nel's method, that is, using 



s 2 \ 



Choosing m = 0, in the first instance, gives the 
Fresnel integral 



which is the same as the largest term in the ex- 
pansion of the cylindrical Bessel function. It 
may be verified that this is true for all values 
of m. 



This demonstrates (20). Substituting into the 



equation (|3| for an eddy gives 



Ae 



-i(ui t+m8—k s s) 



jm(k r (T ) 



k^ds 



(21) 



The substitution is valid in the case k s = 0, and 
the equivalence of the two expressions for k s ^ 
follows by applying a Lorentz boost. 

When m = or 1 , which are the cases of inter- 
est in the article, the Legendre polynomial P m 
reduces to P m (cos9) = cos(m0) and the inte- 
grand is a solution to the wave equation. For 
larger values of m an expansion in terms of Leg- 
endre polynomials is required. 

Given that the integrand is a solution to the 
wave equation, the path of integration need not 
be restricted to a straight line. Any path of in- 
tegration, corresponding to any curving of the 
eddy, can be taken, and the result will still be 
a solution to the wave equation (strictly speak- 
ing, an approximation is involved if k s ^ be- 
cause we have used a Lorentz boost in the above 
derivation, which is restricted to a straight line, 
but we will neglect the perturbation). 

In the coordinates of Figure [3j replacing k s s in 
(21) by n4> because the motion must be single- 



valued gives 



A e 



-iuj Q t 



Rn 



R mn « / e^ me '- n ^ ]m (k r a) R c 



which are the equations used in the article. 



j m (cr) ds ~ - / sin ( r + — ] ds (r — > oo) 



s 2 



D Spin symmetry 



which evaluates to 



The main text shows that the density pattern 

associated with a Rn sonon has the symmetry 

of the Pauli matrices. A possible physical inter- 

/~2~ / tt\ pretation of the associated spin-half symmetry 

Jm (a)ds^-J— cos(r + -J (r ^ oo) is as fo u ows . 
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It is shown in the main text that two iden- 
tical sonons have a preferred alignment where 
there is an antinode of density between them 
and their spins are opposed. Refer to Figure 
[8j which considers the effect of rotating such a 
pair in space whilst maintaining this alignment. 
It may be noted that one complete rotation as 
drawn reverses the sign of the density pattern. 
This symmetry may be relatred to the spin-half 
symmetry associated with the Pauli matrices. 




Figure 8: Two Rn sonons maintaining their 
preferred alignment (with an antinode between 
them), drawn as if they were close to one an- 
other for clarity. A complete rotation is pro- 
duced by the following idealised steps: (a) the 
blue sonon remains stationary whilst the white 
one rolls around it by half a turn, as if locked 
by gears, producing (b). The entire structure 
then rotates rigidly through a further half turn, 
producing (c). Observe the sign is reversed after 
what is effectively a complete rotation. 



E Dirac equation 

The expressions related to the Pauli matrices in 
the main text can be written in Cartesian coor- 
dinates at large r 

3$ 



% a 



j ■ 



dxj 



±JL$ 



(22) 



where <x,- are the Pauli matrices and summation 
over j is implied. The contributions from the an- 



gular components have been neglected because 
of the term - which appears when transforming 
the angular components from spherical to Carte- 
sian coordinates. 



Notice that (22) is analogous to the Dirac 



equation in the case of a stationary particle, 
where the time derivative is assumed to vanish. 

The main text shows that the irrotational so- 
lutions to Euler's equation are Lorentz covariant 
to a very good approximation at all amplitudes. 



Consequently (22 ) can be generalised to an equa- 



tion which applies at all velocities, by writing it 
in a Lorentz covariant form (using the procedure 
which is shown in detail in the main text for a 
scalar field). It may be noted that there is an 
analogue with the Dirac equation, which is also 



Lorentz covariant and reduces to (22) when sta- 
tionary. 

F Fluid energy near a sonon 

Here we obtain an approximate expression for 
the kinetic and potential energy in the fluid mo- 
tion within a sphere of radius R of a sonon. 

In the main text it is shown that the density 
pattern at large distance from a sonon is given 
approximately by 



R,, 



R B.$ 



where B is of order 1. Here we will calculate the 
energy for the term $ = sin (k r r)/r (the term 
in cos(k r r)/r behaves the same), replacing the 
constants which are of order 1, A and B, simply 
by 1. 

The amplitude of these standing waves is 
Ap = p R /r at their peak. Substituting into 
the fluid speed [ID] , which is cAp/p Q , gives 



Upeak 



R 



The kinetic energy density p u 2 /2 is equalled 
by the potential energy density when averaged 
over a cycle, and so the total energy inside a 
sphere of radius R is 



E 



fluid 



4irr 2 p u 2 dr 
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The spatial integral gives a factor of 1/2, and 
averaging over time gives another factor of 1/2 
giving 

Efiuid ~ vRl R PoC 2 
which is the equation used in the main text. 



G Boundary condition: 
coupling mechanism 

In the main text, it is shown that the energy in 
the waves surrounding a sonon is given by ([9]), 
which will be minimised if sonons align them- 
selves so their waves interfere destructively at 
large distance. Here we describe a coupling 
mechanism which allows this lowest energy state 
to be reached. Huygens discovered a similar 
mechanism in pendulum clocks in 1665 [TTj . 

Two similar sonons will behave as oscillators 
which are weakly coupled through he nonlinear 
term in Euler's equation. They will consequently 
have two modes of coupled oscillation, in-phase 
and antiphase, of frequencies 



./ 



fo±" 



(23) 



where v depends on the strength of the coupling. 



If both modes in (23) are present, there will 



be a component of motion at the difference fre- 
quency 2u, which will be coupled to longitudinal 
waves by the nonlinear term in Euler's equa- 
tion. The lower energy state can be reached 
by these waves carrying energy away, which 
will be favoured because of entropy considera- 
tions (assuming a low temperature at large dis- 
tance). The mechanism will be efficient because 
the waves have a frequency which is significantly 
different from the natural frequency f a of the 
nearby sonons of a given type, and so they will 
have a long mean free path. 

Notice that this mechanism will also have the 
effect of locking together the frequencies of the 
sonons of similar types in a local region, even 
if they initially have slightly different sizes and 
frequencies. 



H Redshift near a massive 
body 

In the main text it is shown that the sonons 
near a massive body become enlarged and their 
frequencies redshifted. These are still solutions 
to Euler's equation, since if £(x, t) is a solution 
then so is £(ax, at) where a is a scale factor. In 
units where a = 1 at large distance, then a < 1 
near the body. 

Here we outline a possible relationship be- 
tween this redshift and the gravitational forces 
which are studied in analogue gravity. The treat- 
ment here uses the analogue with quantum me- 
chanical and classical motion which is set out in 
the main text, and it is only intended to apply 
at large distance. 

In a uniform situation, the frequency / of the 
sonons of a particular type at any given position 
will be the average of the frequencies in the vicin- 
ity, due to transport processes, giving V 2 / = 0. 
This applies to all types of sonons and therefore 



V 2 a 







The spherically symmetric solution to this 
equation can be written in the form 

, 2r 

a = 1 

r 

where r Q is to be found. 

A test particle will be perturbed by this field, 
since its frequency will be redshifted by a due to 
the Huygens correlation discussed in £jGj This 
corresponds to the relationship 



to 



u 1 - 



2r D 



and multiplying by h gives an equivalent ex- 
pression for the energy perturbation, which is 
propotional to the natural frequency or inertial 
mass of the sonon and inversely proportional to 
distance. 

At large distance, this perturbation corre- 
sponds to an inverse square force proportional 
to the inertial mass of the test particle. From 
the conservation of momentum there must be an 
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equal force on the attracting object, and, from 
symmetry, it must be proportional to the prod- 
uct of the inertial masses. This is analogous to 
Newton's law of gravitation. The effects of the 
redshift have been neglected and consequently 
this is only true at large distance. 

I Analogue of the Coulomb 
force 

Here we calculate the fluid dynamical force be- 
tween two stationary i? n sonons when they are 
in their preferred alignment shown in Figure [5} 
We will also use the calculation of the fluid 
dynamical force between two identical sources 
or sinks of magnitude Q which are distance r 
apart, which is found in Faber's textbook [HH 
§4.5]:- 



Po 



Q 2 



° Arrr 2 

The flow pattern near a Rn sonon is sketched in 
Figure |4^a). There is a region where the flows 
reinforce and the nonlinear ity is significant. This 
will be modelled el. S 9 cylinder as shown in Figure 
|4]^b). Density variations are rolled up into the 
effective mean flow speed. 

In this idealisation, the external flow derives 
from a source and a sink at the ends of the cylin- 
der in Figure |4](b) , each of magnitude 



Qi 



±7T 



aR a \ (5c 



(24) 



These sources and sinks rotate with the sonon 
at angular speed u D . When viewed from large 
distance in the x direction the flows cancel one 
another at the instant shown, as may be seen 
from the symmetry in Figure |4} However, when 
viewed from large positive y there will be incom- 
plete cancellation at the instant shown because 
of the time delays associated with information 
crossing the structure as it rotates. The incom- 
plete cancellation is quantified by the approxi- 
mate factor 



COS bJ D 



t 



bR n 



— COS U). 



t + 



bR n 



and the resulting effective source becomes 
7ra 2 R 2 

Q = ° (3c 2 sm(u t) sm(uj bR /c) (25) 



Substituting into (14) gives the force between 
these sonons 



Po 

64r 2 



! r > 1 <2 2 • 2 , M 2 / ^ WoRo 

na RJ3 c sin (uj t) sin 



Now, sin 2 (o; t) averages to 1/2 over a cycle, and 
the volume of the cylinder is V — ira 2 bR? /2, 
giving the mean force 

r '.' p n V0 2 c 2 R n S in 2 ( b -^-) (26) 



1 <r 
64r 2 b 

The force produces an acceleration a = F/m 
where m is the mass of the fluid in the cylin- 
der (where the nonlinearities are concentrated), 
which can be estimated as p Q V neglecting den- 
sity fluctuations. In the corresponding calcula- 
tion of the attraction between bubbles which are 
subjected to ultrasonic vibration, the motion of 
the displaced fluid and the associated dressed 
mass is also taken into account [10J. This would 
give an additional factor of order 1/2, but we 
neglect it for simplicity, and it is likely to com- 
pensate for our neglecting density fluctuations. 
Thus we will use the approximate value 

F 

PoV 

conventional units, F r = 



in 



Expressed 
(hu/c 2 )a. Substituting into (26) gives 

a 2 (3 2 fthj n R n . r, ( b u n R, 



646 



sin 



(27) 



For a linear eddy, u 2 = c 2 (k 2 + k 2 ), as dis- 
cussed in the main text. For a Rn sonon, which 
has one twist around the torus, k s = l/R ap- 
proximately, and so 



U) n R n 



Vk 2 R 2 + 1 



Substituting into (27) gives 

a 2 (3 2 



C sin 2 (C) 



he 



646 

where C = \Jk 2 R 2 + 1. This is used in the main 
text, dropping the unnecessary suffix c. 
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J Analogue of Maxwell's 
equations 



The inverse square force in (15) has been cal- 
culated for stationary quasiparticles. Here we 
extend it to moving quasiparticles by using the 
symmetry that the solutions are Lorentz covari- 
ant. Our calculation is analogous to the usual 
extension of the Coulomb force to all veloci- 
ties, and it will be no surprise that the result- 
ing equations are analogous to Maxwell's equa- 
tions. Whilst this function is formal, this section 
may nevertheless be of interest because, along 
the way, it also obtains the Lagrangian from the 
properties of a sonon in a transparent way, and 
it makes explicit the analogue between the phase 
of a sonon and the quantum mechanical phase 
which is measured in a superconductor using the 
Josephson effect. 



The force in (15) is curl-free, so it can be writ- 
ten as a potential gradient, F = gV$. The 
quantity q is analogous to the charge on a test 
particle and $ to the sum of the potentials due 
to all other charges. 

Using Newton's second law in the form 



-Vw 



dk 



gives gV$ = hdk/dt = hV(dS/dt) where S is 
the phase of the sonon, which integrates to 



dS 



T 



constant 



(2f 



This key equation is analogous to the phase evo- 
lution in a charged superconductor, which obeys 
the same equation with q = — 2e and which is 
measured experimentally using the Josephson ef- 
fect 

A sonon moving between two fixed events suf- 
fers a phase change which must be stationary 
with respect to path variations. There would 
be destructive interference if this were not so. 
In free space the phase change is AS = E r/h 
where r is the proper time measured in a co- 
moving frame. In a stationary frame, hAS = 
J (E D + T)dt, where T = (7 — 1)E D is the kinetic 
energy. 



This motion will be perturbed by the interac- 



tion in (28). Defining the Lagrangian 



L = T - g$ 

then HAS = J(E + L)dt. If this phase dif- 
ference is stationary with respect to path vari- 
ations, 5 (HAS) = 0, then the Euler-Lagrange 
relationship d/dt(dL/dvi) = dL/dxi reduces to 



dp 

~dt 



which recovers Newton's second law for a par- 
ticle of charge q near stationary electrostatic 
charges. 



Equation (28) can be extended to a moving 



system by exploiting the symmetry that the so- 
lutions are Lorentz covariant. The Lorentz co- 
variant 4-vector is 



\dS_ 



vs 



Q (® 

h 



+ r 



Notice that A is analogous to the magnetic 
vector potential and Y the gauge. In a coherent 
system, the phase change around a loop must be 
quantised, so that 



VSdl 



Inn 



where n is an integer. This is analogous to the 
flux quantisation which is observed in a super- 
conductor. 



The force in (15) is inverse square, and in gen- 



eral such forces obey, in suitable units, 

-e D V 2 $ = Pc 

where the charge density p c represents the 
sources. This equation can also be extended into 
a Lorentz covariant form in the usual way, giving 

.0 1 d 2 



c 2 df 



Pc 

e 



A 



Likewise, extending the density to a 4-current 
gives, more generally, 

c 2 dt 2 

where J is the current. These are analogous 
to Maxwell's equations for the electromagnetic 
field. 



p 3 



15 



K Radiation from an oscil- 
lating sonon 

As discussed in the main text, the density pat- 
tern of a sonon at large distance is given, without 
loss of generality, by 

f = ABi e- tWot — sm(k r r) (29) 
r 

If a sonon oscillates so its position is x = 
x sm(Qt) where x Q is small and its velocity of 
motion is much less than c, then this becomes 



f = AB 1 e- tuJot — sin[A; r (r + r')l 
r 



where 



r' = x o cos(0) sin Q (t — - j 



and is the angle from the direction of oscilla- 
tion. This may be written in the form £ = Xwipw 
where for low amplitude motion, using the first 
term in the Taylor series expansion in k r x 



I AR ° I- 



B x e~ iWot cos(k r r) cos0 



Notice that we have chosen to distribute the fac- 
tors so that ip w is a solution to the wave equation 
in its own right. This distribution differs (on a 
technicality) from that used for a sonon. These 
are the equations used in the main text. 



and uj /D respectively, where the Doppler shift 
is 



D = 



c + v 



V 



The total energy in this case is hoj (D + 
l/D)/2 = ^ftiOo, which matches the energy of 
the moving particle. A similar calculation for 
perpendicular motion shows this equivalence to 
be independent of direction. 

A similar consideration applies to momentum. 
Suppose the momentum of each wave is p w = 
(5E W = /3hu Q where f3 is to be determined. When 
at rest, the waves have opposite momenta so 
the total momentum vanishes. A Lorentz boost 
changes this to p — (3hu (D — l/D) = 2/3Ev. 
This must be the momentum of the particle, 
p = Ev/c 2 , and therefore 2/3 = 1/c 2 and 



Pv 



E„ 



L Energy and momentum 
of a wave 

This appendix calculates the ratio of energy to 
momentum for a wave associated with a particle. 

When at rest, the energy of a particle is E — 
hu so each of the outgoing and incoming waves 
must have energy E w = hu /2. This equivalence 
also applies after a Lorentz boost. In the case of 
parallel motion, the frequencies of the outgoing 
and incoming waves are Doppler shifted to u Q D 
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